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Periodic thick-brane configurations and their stability 
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I investigate models with scalar fields in 5 dimensions that exhibit thick-brane configurations 
with a non-trivial metric. I show that an appropriate coupling to the scalar curvature allows for 
periodic configurations, which, however, are unstable under small harmonic perturbations. A model 
for stabilizing these configurations is proposed and discussed. 
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^H ■ I. INTRODUCTION 

o: 

^ . One of the most interesting proposed extensions of the standard model is based on a 5-dimensional space-time [1] 

^ ' containing one or more branes, 4-dimensional subspaces where some of the fields are localized. Such models usually 

i— 5 I assume that space time has a topology M x B, where M denotes the usual Minkowski space and IB is a manifold that 

ly-N . may or may not be compact, and a non-separable metric (I consider here only the so-called "small extra dimension" 

models); the original models assumed IB to be one-dimensional but this has since been generalized [3]. When this 

' ^ ', space-time configuration and brane content correspond to a stable solution of the field equations, it can serve as a 

"^j ' vacuum where quantum fluctuations propagate; the assumption being that this background configuration corresponds 

JL , to a minimum of the effective action for the complete system. 

Qj • This paradigm offers an innovative solution to the hierarchy problem and a wealth of new effects that may be 

(~| [ observed at the LHC. Motivated by this, several authors have constructed realistic or semi-realistic models based on 
this idea, and have studied a variety experimental signatures [2] . 

The stability of the background configurations under small harmonic perturbations has also been studied [1, 7, 
^ . 25]. In the seminal publications [1] it was shown that for the models there considered most perturbations are not 
■»H- ' destabilizing; the one exception being the dilaton, a self-similar scaling perturbation which is neutral (has no quadratic 
^^ , potential term or, equivalently, is associated with zero frequency). This problem, however can be eliminated by the 
0^ • addition of appropriate bulk scalar fields [8], which has been studied in a variety of cases: for single scalar fields 
^^ ' coupled minimally to gravity [8, 9], non- minimally coupled scalars [10], and for Brans-Dicke theories [11]. 
f~^ ' These higher-dimensional theories, being non-renormalizable, have an intrinsic ultraviolet cutoff scale A beyond 

^^ I which they are not reliable (at least within perturbation theory). In models where branes are infinitely thin (in the 
'""' . fifth coordinate), it is tacitly assumed that their structure will become manifest at scales <; A. But this need not be 
the ease, and a certain amount of attention has been paid to the possibility that the dynamics at scales below A is 
responsible for the brane configuration. A variety of such models have been studied in the literature for non-compact 
IB with one [12-14], or several [15, 16] minimally-coupled scalars. In addition, the stability of such models has been 
^ , investigated for the case of non-compact B and one [13, 14] or several [16] minimally coupled scalars, and for B 
j^ ■ compact and a flat metric [18]. 

Models with B compact that exhibit periodic, stable, brane-like configurations generated by physics below A 
("thick branes"), and which have a non-trivial metric are more difficult to construct. The reasons are, first, that for 
the simplest models the background configurations cannot satisfy the periodicity requirement (e.g. they fail to satisfy 
one or more of the sum rules listed in [17], as required for consistency). And second, general considerations [18, 19] 
apparently preclude the stabilization of thick-brane solutions using only scalar fields. 

The goal of the present paper is to provide mechanisms that overcome these two difficulties. Specifically, to 
exhibit a class of 5-dimensional models containing gravity, scalars, and antisymmetric tensor fields, which admit 
periodic kink-like solutions that are perturbatively stable at scales below A. Though the theories thus obtained have 
no phenomenological applications {e.g. they exhibit but a small amount of warping, so that no significant mass 
hierarchy can be generated), they are of interest because they can address these two problems. Though there are 
some indications that more realistic configurations can also be obtained (sect. II A), I will not attempt to construct a 
phenomenologically viable model, nor will I not attempt to address the nuich more ambitious and difficult problems 
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of confining fields to the branes and of global stability. For a different class of theories, based on the "large extra 
dimensions" paradigm [5], there are realistic models that include both dynamically generated branes and a confining 
mechanism for the Standard Model fields; see for example [6]. 

The calculations presented are essentially classical, amounting to obtaining solutions to the equations of motion 
that are stable under harmonic perturbations. The usefulness of these results lie in the well-established connection 
between such solutions and related quantum objects [20]; in particular, the classical background solutions are to 
be interpreted as the lowest order semi-classical approximation to a quantum vacuum, and the frequencies of the 
harmonic perturbations as the energies of the low-lying excitations. 

The plan of this paper is the following. The next section introduces a class of models involving only scalars, 
that support periodic configurations, which in favorable cases are similar to the ones in [1]. These background 
configurations are unstable (Sect. Ill), but this problem can be solved through the introduction of antisymmetric 
tensor fields (Sect. IV), adequately coupled. The last section contains some parting comments and observations, 
while some mathematical considerations are delegated to the appendices. 

II. THE MODEL (BASIC VERSION) 

I first consider the problem of constructing kink-like configurations that involve gravity and are periodic in the 
fifth coordinate of a 5-dimensional space-time^ . Though it is well known [22] that scalar models with Lagrangians of 
the type £ ^ (90)^ — V{(j)) do allow stable kink and multi-kink configurations, none of these satisfy the periodicity 
constraint (the simplest way of seeing this is by noting that such configurations violate one or more of the sum rules of 
Rcf [17], as shown in that same publication); fortunately a simple and natural modification overcomes this obstacle. 
I consider models of the type [10], 

C = 2M|ii? - -g'^4>r,i4>r,j - ^^rs{4>r + flr^){k + ns^)R + V ((j)^ , (1) 

where Afpi denotes the 5-dimensional Planck mass, and R the scalar curvature, a comma denotes an ordinary deriva- 
tive, is an iV component real scalar field with components 0,., and ^ is a real, symmetric N x N matrix. The vector 
n denotes a constant direction in field space, while 1^9 is a scale of the same order as Mp^ . The sign of the potential 
term is chosen for later convenience. 

Models of the type (1) have been studied in the literature, in particular Ref. [18] provides fiat-space, single-field 
examples where the background solutions are both periodic and stable; unfortunately these models do not extend 
easily to the case of a non-trivial metric. 

I assume now that Mpi is the largest scale in the theory and consider solutions that have a small but non-trivial 
deviation from a fiat metric (recalling that ip ^ Afp( ): 

= (j) + 0{l/^); 
5y = Vij + -hj + -^hj + 0(l/i/) ; (2) 

ip Lp'^ 

where rjij = diag(— 1, 1, 1, 1, 1) denotes the flat-space metric. Substituting in (1) gives, after some algebra. 
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g^ ^/i'^''%,fe - h^\kh..,'^ ~ p^q' ] + 0{l/p) , (3) 



where space-time indices are raised and lowered using the fiat metric 77, and where 



In the following i,j, etc denote 5-dimensional space-time indices, Greek indices will refer to the 4-dimensional non-compact coordinates; 
the metric signature is(— 1,-|-1, +!,+!, -(-1). I use the conventions of Landau and Lifshitz [21] for the definition of the Riemann and 
associated tensors. I denote the compact coordinate by y, and the non-compact ones by x^, ^ = 0, 1, 2, 3. The indices r, s, etc. label 
the scalar-field components. 
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The corresponding equations of motion are: 
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which can also be obtained by expanding the Einstein and field equations and using (2). 



A. Background solutions 



With the goal of preserving 4-dimensional Lorentz invariance I look for background configurations of the form 



4> = *(y), 

hij = Hij{y), H^y ^ -2a{y)rii_,^, H. 
kij = Ki]{y), K^ 



4m 



fj.iy 
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0; 
0; 



(6) 



which solve (5) to lowest order in Lp provided 
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where a prime denotes a y derivative. 

The last two equations in (7) describe the zero-energy, classical, non-relativistic motion in a potential V; the mass (cf. 
eq. 4) equals one except when $ is parallel to a, in which case the mass is 1 — (16/3)^ (for a different interpretation 
when ^ > 3/16, sec A). Higher-order corrections (in l/(p) can be determined similarly, but the expressions are 
cumbersome; an example is presented in B. 

The warp factor in this type of models is ~ 1 — 2(j/lp, and cannot generate a significant mass hierarchy. Note however 
that the expression in B for r shows that this function is also periodic, with the same period as cr; which suggests 
that periodic solutions exist also for moderate values of Lp. Verification of this conjecture, as well as a determination 
of the range in Lp for which it holds, lies beyond the scope of this paper. 

In the following I will choose coordinates in field space such that 



5rs => m = diag ( 1- —^, 1,1,1, • 



(8) 



"^1 = 1 - yt, 



and will denote by 771^ the eigenvalues of 777: 

1 a 

rrir = 1, 2<r <N (9) 

A judicious choice of the potential V and initial conditions will lead to solutions that are periodic in y, 

*(7/ + L) = *(7/), (10) 

with L determined by V and the initial conditions; a will also be periodic with the same period. For example, adopting 
the basis (8) and taking 



V{4>) = ^mrUri(t>r) , 



(11) 



the equations (7) yield 



a = 2$i/(3ni) ; $/' = /,($,); (12) 

where fr = —dur{4'r)l d4>r (no sum over r). An appropriate choice of Ur and initial conditions will then generate 
solutions obeying (10). 

Adopting such solutions as background configurations allows the identification y = y mod L, which amounts to 
a compactification of the fifth dimension. Such configurations can serve as vacua for the theory, provided they are 
stable. 



Sum rules 



It is of interest to see how the present models avoid the obstacles listed in Ref. [17]. The most severe constraint is 
obtained by assuming the existence of periodic solutions and integrating the second equation in (7) over a period: 



® dy^'m^' = 0. 



(13) 



When ^ = this can be satisfied only when $' = (in the presence of localized branes there are additional contri- 
butions and the sum-rule can be satisfied by non-trivial $ configurations [17]). In contrast, when ^ 7^ 0, non-trivial 
solutions are allowed provided m is not positive definite; from (8) this corresponds to 



F 3 
^>16' 

which I assume henceforth. In this case the mass matrix (8) takes the form 

m=|m|S'; S ^ dia.g{-l,+l,+l,- ■ ■) ■ 



(14) 



(15) 



C. Comparison with the 2-brane Randall-Sundrum model 

In the 2-brane Randall-Sundrum model [1] the background of the form (6) is also adopted, but instead of introducing 
scalar fields it is assumed that the space-time contains two branes with cosmological constants ±A, and a bulk 
cosmological constant A < 0. The resulting field equations are 



24Af|icr'2 = -Aip^ ., 



12i\/|i(T" = Xip [S{y) - S{y - L/2)] , 



(16) 



where < y < L is the range of the compact coordinate. 

In order to compare this to the previous results I again adopt (8) and assume the potential takes the form (11) 
Then, from (7), 
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Using (12), it is clear that for a very fiat potential ui (see Fig. 1), ^'^ will be almost constant except at the turning 
points where it will rapidly drop to zero, while $" will be almost zero except at the turning points, where it will be 
large. This type of potential then yields configurations qualitatively similar to those derived from (16) for the case 
A<0. 
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FIG. 1: Scalar potential leading to configurations similar to those described in [1] 



III. STABILITY 

The usefulness of the above solutions as background configurations depends on their stability; the minimal require- 
ment being that all periodic solutions to the linearized perturbation equations are bounded in time, up to coordinate 
transformations. As usual, such linear perturbations can be assumed to depend harmonically on the non-compact 
coordinates. 

The relevance of coordinate transformations can be illustrated by considering a background solution of period L, 

*(2/) = J2 f^e^™'^/^ ; 'yiy) = E ^^e^^™^/^ ; (18) 

and a perturbation that consists in the replacement L — > (1 — e)L: 

S^cjy = [*]i^i_,i - * = 6J/*' -eLY, e^-^^/^a^f^ , 

hgi^L = (5^.944 = . (19) 

Though this suggests the need to include periodic perturbations multiplied by linear functions of y, this is not the 
case since these modes can be absorbed using appropriate coordinate transformations. Under x* — > a;* 4- dx)j 

Scoordgij = Ci:j + ij;i+ 0{£, ), 
'JcoorcK/' = e^^^ + 0{e), (20) 

where a semicolon denotes a covariant derivative, and choosing ^4 = —ey, C^^ = 0, a combination of (19) and (20) 
gives 

S'l944 = -e ; (5^54^ = ; 

S'^g^u = -2eL^e2""^/^9LS„77^,, 
6'^<f> = -6i^e^™^/^aif„, (21) 

where 6']^ = Sl + Scooid- The variations 5'^gij, S'j^cp BJ:e then equivalent to (19) and arc periodic in y. 
Hence I will look for solutions to (5) of the form 

= <^{y)+e'P-''\m\-^/^x{y), p■x^p^x^'■ 
K, = i7„(2/) + e*P-7,,(y), (22) 

keeping only first-order terms in the perturbations x, 7ij (the matrix |to|^^" - cf. eq. 15 - is included for later 
convenience) . 

In order to simplify the calculations it is convenient to choose coordinates such that 

74^ = 0; 7l=0. (23) 

(as before, indices are raised and lowered using the flat-space metric). These conditions, however, do not completely 
fix the coordinates: using (20) with 

e = iX'{y)e'P-^, e = X{y)p>^e^P-^ ■ X"+p'X = 0, X^l/^p, (24) 

yields 7^11/ — > 7^^ + 2ipf^p^X and 744 — > 744 + 2iX" , which preserve (23). 
Expanding 

^A^- = ^^i + (^^, - ^^ Tt + %.; p^p^'l^.u = ri'"'%. = , (25) 



where p'^ = ^l^'^P^iPu, and substituting (22) in (5), yields 

TT = -'d, Vl+ 744 = i-d , 

%. - P^%.^ = , 744 + P'744 = M" - pH , 

and 

d^ 



(26) 



i?oX = Vx; Ho = -j^ + SU, (27) 



where the first equation in (26) enforces the second condition in (23), S is defined in (15), and 

4£ 1 / d'^V \ 

3|a| ^/\mrms\ \o4>r Oc/Js J ^^^ 

Destabihzing modes correspond to periodic (in y) solutions to these equations when p^ > 0; in this case (26) require 
7^1/ = 0; as otherwise the solutions would not be periodic. Note also that the modes x ~ ItJ.^ — Tr = 0, T l = —744, 
with 744 obeying the homogeneous equation 744 + ^^744 = 0, are apparently destabilizing, but these modes are 
coordinate artifacts that can be eliminated by a transformation of the form (24). 

Since the homogeneous solution to the 744 equation in (26) can be eliminated using (24), we can assume that 744 
is determined by 1? oc a • %. The stability of the background configuration then depends only on whether (27), the 
equation satisfied by x, has periodic solutions only for p^ < Q. The analysis of this equation is complicated by the 
fact that SU is not Hcrmitian (for the usual definition of the inner product), so that its similarity to the Schrodinger 
equation is not very useful in this case ^. Instead I argue as follows 

The general solution to the equation m$" = F in (7) depends on 2N constants {e.g. the values of $ and *' 
at y = 0). The periodicity requirement (10) introduces N restrictions and one new parameter, since L is not fixed 
a priori. Finally, (13) imposes one additional constraint. As a result the general periodic solution will depend on 
N arbitrary parameters (including L); since the equation of motion does not depends explicitly on y, one of these 
parameters can be taken as some initial value j/o- The general periodic solution can then be written 

* = *(?/- 2/0; c;L), c= (ci,...,cw-2)- (29) 

Let now 

X^°)=a,„$, x^^-i) = i^^*, x''^=5e.*, i = l,...,7V-2; (30) 

where 5'j^ is defined below (21). Then, by taking the appropriate derivatives of the equation of motion (7), it is easy 
to see that 

iJoxW=o, i = Q,...N-l. (31) 

Since the constants j/Oi L and c are independent, the modes x i i = 0, . . . , N ^ 1 will be linearly independent (and 
periodic). 

The equation HqX = is a linear, second-order differential equation for the component functions Xr, so the general 
solution will contain 2N arbitrary constants. The periodicity condition, however, imposes N constraints (note that 
now L is fixed by the background solution), so that the general periodic solution will contain only N independent 
constants. Then, since the equation is linear, any periodic solution to (27) with p^ = can be written as a linear 
combination of the x'*'. 

Consider next those background solutions that remain near an extremum $cxtr of the potential V (the presence of 
an extremum is a necessary condition for the existence of periodic solutions (10); when mi < 0, as for the cases of 
interest, this extremum is a saddle point). For these cases the relevant equations are, approximately, 

*" = SU^ , -x" + SUx = -p"-x ; (32) 

where 

!>, = |m,|i/2(*-*,,t,)^ ; tJrs = -\mrm,\-^^^ ( ^-^] . (33) 



'^rVs J A.— 



Now, for there to be periodic solutions $, SU must have at least one negative eigenvalue — w^ (w real); denoting by 
V the corresponding eigenvector, it follows that the y-independent choice x = v is a solution to the x equation with 
p^ = w^ > 0, which is then a destabilizing mode. It follows that all small amplitude solutions are unstable. 

Now consider a small amplitude solution $(j/ — yg;c';L'). another arbitrary solution $(?; — ?/q;c";L"), and a 
smooth path in parameter space {yo{s),c{s), L{s)}, < s < 1, such that {yo(0), c(0), L(0)} ~ {yQ,c',L'} and 



^ Hq is Herinitian under the inner product (xi|X2)5 = f dyXiSx2t but, thought this can be used to show that p^ must be real, it does 
not provide information about its sign since this inner product is not positive definite, and the arguments based on the Sturm-Liouville 
theory [18, 23] are not applicable. 



{yo{l),c{l),L{l)} = {j/g,c",L"}. Each value of s defines a periodic solution $ with period L{s), form which a 
corresponding U is constructed. As a result all the modes % will depend smoothly on s, and so will the eigenvalues 
—p^; in particular, the p^ = subspace is N dimensional for all s. For s = we know there is at least one mode with 
p^ > 0, but then this mode must remain destabilizing for all s, for if it were to change from destabilizing to stable, 
there would be a value s at which its eigenvalue p^ would vanish, so that for s = s the p^ = subspace would have 
dimension A^ + 1, which we saw above is impossible. It follows that the corresponding solution $(?/ — yQ;c";L") is 
also unstable, and since the parameters {yQ,c" , L"} are arbitrary, it follows that all background solutions in the pure 
scalar-gravity model are unstable. 

There is a subtlety in the above argument. One can easily imagine potentials with two or more extrema, each with 
different number of destabilizing modes, yet the above argument seems to indicate that the number of destabilizing 
modes cannot change. This apparent contradiction is resolved by noting that there is no periodic solution that can 
interpolate between the corresponding small-amplitude solutions. Imagine a potential with 2 extrema and small- 
amplitude solution near each, denoted by $^^^(j/) and $'^^^(y), and for which the corresponding x equations in (32) 
have different numbers of destabilizing modes. Then for any periodic function $(y, s), < s < 1, of period L[s), that 
solves (7) for each s, and such that #'^^^ = #|s=o, ^^^"^ = ^|s=i, there will be an intermediate value < Soo < 1 such 
that L{s) — > oo as s — > Soc- Thus, to every extremum one can associated a "region of influence" determined by all 
solutions that can be reached by an interpolation with L{s) finite for all s; the above argument shows that all such 
solutions are unstable. I will assume that all periodic solutions can be characterized in this way, that is, that any 
periodic solution can be "deformed" into a small amplitude solution near an extremum of V while keeping the period 
finite. These results extend the arguments of [18, 19] to the class of models considered here. 

Similar considerations can be followed in case the potential is not quadratic near the extremum. 

IV. ANTISYMMETRIC TENSOR STABILIZATION MECHANISM 

The instability of the solutions to (7) is reminiscent of the well-known instability of soliton-like solutions in more 
than 1-1-1 dimensions [22]. In the soliton case stable solutions are obtained by introducing gauge fields, here I will 
pursue a different approach based on the introduction of an antisymmetric tensor field A that can propagate in the 
bulk, and which has the following interactions: 

r, — n^J n'^^ n"^"^ A , 4, — - k^ n^J n'^'- A -, A , 4- - Kn^^ n'^'' T, A , ■ 

'-A — TT.y 9 g -^ikm^jln ^1^ 9 9 ^ik^jl + ^1^9 9 'Jik^jl ; 

Jik = 2J l"^"""^*! ''^''*'^'">*'?^«>'= ' (^^) 

where Aijk = 5Auj^ and Xrs + Xsr ~ (the couplings A could depend on also, but such a case will not be 
considered here). When considering the quantum aspects of this model it is convenient to rewrite Ca using an 
auxiliary Stiickelberg-like vector field (see, for example, [26]); the expression (34) then corresponds to the "unitary" 
gauge where this auxiliary field vanishes. The presence of higher-derivative terms, while innocuous classically, can 
have dire quantum effects; still, since these theories have an intrinsic UV cutoff, these are avoided with appropriate 
constraints on the scales related to the couplings A (see sect. V). Antisymmetric tensors have been considered 
previously in the literature in RS-like models with or without fundamental scalars [27], in non-periodic thick-brane 
models with weak scalar-tensor couplings [28], and in exotic-Lagrangian models addressing the self-tuning of the 
cosmological constant in RS-like models [29] . The type of models described by (34) has apparently not been discussed 
previously within the present context. 

I will assume that k = 0{ip'^^^) so that A ^ 1/(^9; then 

Ca = -\{kA- j),,{kA - jy^ + ^j,,r^ + o{i/^) , (35) 

where indices are raised and lowered using the fiat space metric r]ij . In this case the field A can be integrated out and 
one can work instead with the effective Lagrangian 

CeS = \AjJ'' , (36) 

up to corrections of order 1/ip. 

The addition of C^s does not change the background equations (to lowest order) so that the results of sections 
IIA-IIC are not modified; in terms of the antisymmetric tensor field this implies that it vanishes in the background 
configuration. 



A. Stability with antisymmetric tensors 

The perturbation equation is modified by the addition of (36). Instead of (27) one gets 

HoX = -p'[x-SB{B-x}], B. =^|m,|i/2A,,$;. (37) 

s 

As before the perturbative stabihty of the background configuration $ is guaranteed if this equation has no solutions 
for p2 > 0. 

Given the form of (37), it is clear that the modes x'*' defined in (30) again provide N solutions corresponding 
to p^ = 0. The argument of section III then implies that the background configurations will be stable provided the 
small-amplitude background configurations near an extremum arc stable; this can be realized when B 7^ 0. 

For background configurations that remain near an extremum of the potential V we still have ^" = Stj^ as in 
(32), but the equation for x becomes 

-x" + SUx = -p'x+p'SB{B-x), B.-5]A,,$',. (38) 

s 

In order to have periodic solutions $, SU must have one or more negative eigenvalues — cj^, with v^ the corresponding 
eigenvectors; however, the constant modes X = Va are not solutions of (38) when B ^ 0. 

I will now consider the case where the potential takes the form (11), with Ur of the shape given in Fig. 1 for all r. 
Then each $^ will be approximately a square wave, so that I can approximate {^'r)^ — > a^ = constant, and the sum 
rule (13) becomes 



a-i 



^a^. (39) 



r>l 



In addition, {SU)rs = Srs^r is diagonal, and each V,. can be approximated by a delta-function comb, with negative 
magnitude and half the period of $r- 

I will now make the simplifying assumption that the periods of $r for r > 1 arc much smaller than that of $1; for 
example, if i denotes the period of $1, that of $r, r > 1 can be chosen to be l/rir with rir a large prime number; the 
period L in (10) then equals (J| rir) £. In this case SU and B (g) B in (38) will have entries that are slowly varying and 
others that vary very rapidly, and this equation can be treated using the procedure described in [30] : up to corrections 
of order l/rir, the solution can be obtained by considering only the slowly varying terms in (38): 

-Xr - (^r,! VlXl = -p'^Xr + P^ Sr ^ BrsXs (40) 

s 

(no sum over r) where, as before. Si = —1, S^yi = -1-1, and 



Brs — 2_^ (-^rlAsl + Xru^su) du ■ (41) 

U>1 

The term containing B in (40) is generated by averaging the rapidly varying terms in (38) over periods <ti (■■ 

The equation for Xr can now by analyzed using standard techniques (see e.g. [31]): denoting X = [XiX'Y there 
will be a non-singular matrix T such that X(j/ + L) ^ Tli.{y), and the system will support destabilizing modes if, 
for some positive p^, T has an eigenvalue equal to one (which corresponds to x being periodic). In determining the 
conditions under which this occurs one is interested only in modes for which |p^| < A^, or, cquivalcntly, [p£ / (Att)]'^ < 
[AL / {AiTAf)]'^ = p, where A/" = H "-r s-i^d i/rir is the period of ^r, defined above; with appropriate choice of these 
parameters one can insure p = 0(10). In this case there are ranges of parameters where not destabilizing modes 
are present; Fig. 2 contains an example for fixed A^s, Or and Vi where the matrix T ^ 1 for p^ > 0, while Fig. 3 
gives, for fixed A^.s, the values of 02,3 where no destabilizing modes occur (within the allowed p^ range). Background 
configurations for this type of models are then stable in such cases. 

Stable solutions occur when the number of scalar fields is > 3; I found no stable configurations in more economical 
models. 

V. COMMENTS 

The main purpose of the present paper was to provide mechanisms that can overcome the periodicity and stability 
problems of brane-like configurations produced by scalar fields. The model defined by (1) and (34) meets these 
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FIG. 2: Value of log |det(T - 1)| when A/" = 1 (left) and TV = 60 (right) In both cases A12 = -0.2, A13 = -0.4, A23 = 0.3, 02 = 
0.05, as = 0.7, and the amplitude of Vi equals 287r/^. Periodic solutions occur for p^ < and do not correspond to destabilizing 
modes. 




FIG. 3: Values of 02,3 for which (40) has no destabilizing modes. Light gray: 
dark gray: p < 20 and the amplitude of Vi equals Att/I. In both cases A12 = - 



p < 30 and the amplitude of Vi equals 29,-k /t, 
0.2, Ai3 = -0.4, A23 = 0.3 



requirements by introducing a coupling of the scalars to the Ricci scalar and to an antisymmetric tensor field that 
can propagate in the bulk. For appropriate choices of scalar potential the background configurations mimic that of 
the RS model but require the presence of > 3 scalar fields with modes of very different though commensurate periods 
(an investigation of whether this hierarchy can be maintained naturally lies beyond the scope of this paper). 

As mentioned in section I these models contain a UV cutoff scale A; when IB is compact they also contain another 
high-energy scale, the compactification radius L. An immediate concern is whether one can naturally assume A ^ 1/i, 
for otherwise the model cannot accurately describe the dynamics associated with the compact directions. This question 
can be investigated using an extension of naive dimensional analysis [4] . For example, a straightforward estimate for 
5-dimensional gauge theories gives A ^ 247r'^/(7|, where (75 is the 5-dimensional gauge coupling constant (which is 
dimensional). The 4-dimensional gauge coupling is 5 '--^ (75/Vi, whence LA ^ 247r^/g^, which indicates that such 
models are reliable at energies above the compactification scale provided Gtt^ ^ g^ /{Air). For compact B the scale A 
also provides a cutoff for the order of the KK excitations that need to be included in the theory; the effects of higher 
order modes are absorbed into a renormalization of the various operator coefficients. 

When the antisymmetric tensor field is introduced, additional constraints appear, such as, for example, those derived 
from the unitarity constraint in (jxj) — > (fxp scattering, which demands (at tree-level) that the scale associated with A 
be below ~ A/(27r); this restriction also insures that the effective Lagrangian (36) does not generate undesirable poles 
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in the scalar propagators. 

In the above models there are no exponentially destabilizing modes below scale A. There are^ however, zero modes 
obtained by taking derivatives with respect to the parameters of the background configuration (sec Eq. 30). Pertur- 
bations along these zero modes then correspond to small deformations of the background configuration parameters, 
under which its period and orbits suffer small changes, but and do not result in an instability. 

Phcnomenologically viable models constructed along these lines must include, in addition, the possibility of orb- 
ifolding. In this case it is most convenient to assume that the potential allows solutions with definite parity under 
y o —y. This can be implemented without additional complications. 

It is worth noting that the present models cannot be stabilized by the 1-loop effective potential [25] . This is because 
in the present case the tree-level potential U supports at least one destabilizing mode (whose amplitude increases 
exponentially with time), not merely a neutral one. Loop corrections are of course present, but they are subdominant. 

Appendix A: Alternative interpretation of the scalar equations of motion 

When t, > 3/16, the presence of a negative sign in S (cf. eq. 15) suggests an alternative description of the equations 
for $ as a geodesic equation, and that for % as a geodesic deviation equation. 

To see this consider an A^ + 1-dimensional space with coordinates $i, . . . , $7v, 9 and metric 

Trs = diag (mi, TO2, ■ • ■ , m^, 1/V^(*)) ■ (Al) 

I will denote by r,s, etc. the indices corresponding to the first N coordinates, then Trs ~ rrirSrs- The geodesic 
equations associated with this metric are 

where p is an affine parameter. 

The second equation can be immediately integrated: d6/dp — cV where c is a constant; substituting into the 
equation for the $r gives 



that reduces to the second-order equation for $ in (7) provided y is identified with cp/\l. 
The zero-energy condition in (7) corresponds to the requirement that this be a null geodesic: 

where a prime denotes a y derivative. It is also clear that the equation for x is equivalent to the geodesic deviation 
equation associated with the metric (Al). 

Appendix B: Background solutions for models with a single scalar field. 

When there is a single scalar field the background solution takes the form 

= $ + e/99H — 

g^, = (1 - 2a/^ - 2r/^2 ^ • • •) 77^, (Bl) 

with ^44 = 1, .94^ = 0. 

For this simple case the field equation is redundant; while the Einstein equations give 

O(^) : 2C$" - 3ct" = , 

0(1) : Vm ~ im<f>'2 - g (2C$' - ia'f - , 

26" - -r" - ma" + 2$'a' - -J>'^ + ($2)" = Q , 
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(B2) 



(I omitted an additional equation to order 1/</j that determines the correction of order 1/^p^ to the metric), where 



e 



e 



e - 4Af3j/^2 



16- 



(B3) 



The 0(1^9) equation together with the periodicity requirement imply 2^$ — 3(t = constant, whence $ must satisfy 

1 



-m$'^ + y($) = 0. 



(B4) 



It now proves convenient to write 



e = $M(y) 



24 6 \ m 



'<''-^ ''*^ + B(„) 



Substituting into (B2) and using (B4) one finds 



a = — $; A' = —^■, 

3 m 



B" = ^A$'' 
3 



•^ 2 y' 3 



32^2 C(l - C) 



0; 



(B5) 



(B6) 



where a prime denotes a derivative with respect to the argument. These equations can be solved by quadratures; in 
particular, 



4> r^ 



fi4>) = -^[Srn' + ^m~0]J dxj d7 



V(7) 
V{X) 



(B7) 



Higher orders can be dealt with similarly. For more than one field the higher-order corrections to (7) cannot, in 
general, be cast into such comparatively simple expressions. 
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